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GALAXY MERGERS AND DARK MATTER HALO MERGERS IN LCDM: MASS, REDSHIFT, AND 

MASS-RATIO DEPENDENCE 

Kyle R. Stewart\ James S. Bullock^ Elizabeth J. Barton^ Risa H. Wechsler^ 

ABSTRACT 

We employ a high-resolution LCDM N-body simulation to present merger rate predictions for dark matter halos and 
investigate how common merger-related observables for galaxies — such as close pair counts, starburst counts, and the 
morphologically disturbed fraction — likely scale with luminosity, stellar mass, merger mass ratio, and redshift from 
z = to z = 4. We investigate both rate at which subhalos first enter the virial radius of a larger halo (the "infall 
rate"), and the rate at which subhalos become destroyed, losing 90% of the mass they had at infall (the "destruction 
rate"). For both merger rate definitions, we provide a simple 'universal' fitting formula that describes our derived 
merger rates for dark matter halos a function of dark halo mass, merger mass ratio, and redshift, and go on to 
predict galaxy merger rates using number density-matching to associate halos with galaxies. For example, we find that 

' the instantaneous (destruction) merger rate of m/M > 0.3 mass ratio events into typical L > f galaxies follows 
the simple relation dN/dt ~ 0.03(1 + /)Gyr~^ (1 + z)^-^. Despite the rapid increase in merger rate with redshift, 
only a small fraction of > 0.4i» high-redshift galaxies 3% at z = 2) should have experienced a major merger 
[m/M > 0.3) in the very recent past {t < 100 Myr). This suggests that short-lived, merger-induced bursts of star 
formation should not contribute significantly to the global star formation rate at early times, in agreement with several 
observational indications. In contrast, a fairly high fraction (~ 20%) of those z = 2 galaxies should have experienced a 
morphologically transformative merger within a virial dynamical time (~ 500 Myr &t z ~ 2). We compare our results 
to observational merger rate estimates from both morphological indicators and pair-fraction based determinations 
between z = — 2 and show that they are consistent with our predictions. However, we emphasize that great care 
must be made in these comparisons because the predicted observables depend very sensitively on galaxy luminosity, 
redshift, overall mass ratio, and uncertain relaxation timescales for merger remnants. We show that the majority of 
bright galaxies at z = 3 should have undergone a major merger (> 0.3) in the previous 700 Myr and conclude that 

\ mergers almost certainly play an important role in delivering baryons and influencing the kinematic properties of 
Lyman Break Galaxies (LBGs). 

I Subject headings: cosmology: theory — dark matter — galaxies: formation — galaxies: halos — 

methods: A^-body simulations 



I 1. INTRODUCTION 

In the current theory of hierarchical structure forma- 
\ tion (LCDM), dark matter halos and the galaxies within 
them are assembled from the continuous accretion of 
smaller objects (Peebles 1982; Blumcnthal et al. 1984; 
■ Davis ct ai. 1985; Wechsler et al. 2002; Fakhouri & Ma 
; 2008; Stewart et al. 2008; Cole ct al. 2008; Neistcin & 
• Dckel 2008; Wetzel et al. 2009). It is well-established 
that galaxy and halo mergers should be more common 
at high redshift (e.g. Governato et al. 1999; Carlberg 
. et al. 2000; Gottlober et al. 2001; Patton et al. 2002; 
Berrier et al. 2006; Lin et al. 2008; Wetzel et al. 2009), 
but the precise evolution is expected to depend on details 
of the mergers considered. Moreover, it is unclear how 
these mergers manifest themselves in the observed prop- 
erties of high-z galaxies and what role they play in setting 
the properties of galaxies in the local universe. Interest- 
ingly, there arc indications that the familiar bimodality 
of galaxies as disks versus spheroids at z = might be re- 
placed by a categorization of disk-like versus merger-like 
at higher redshift (Forster Schreiber et al. 2006; Genzel 
et al. 2006; Law et al. 2007a; Kriek et al. 2008; Mel- 
bourne et al. 2008; Shapiro et al. 2008; Wright et al. 
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2009), although this shift in the dichotomy of galaxy 
morphologies is by no means robust and requires fur- 
ther study. In this paper, we use N-body simulations to 
provide robust predictions and simple fitting functions 
for dark matter halo merger rates and merger fractions 
as a function of redshift, mass, and mass ratio. We use 
our predictions to address two observable consequences 
of galaxy mergers — merger-driven starbursts and mor- 
phological disturbances — and investigate their evolution 
with redshift. 

The tidal interactions inherent in galaxy mergers pro- 
duce concentrations of gas in the remnant centers. For 
major mergers (rn/M > 0.3), models predict that this 
effect results in a significant burst of increased star for- 
mation rate (SFR) compared to the central galaxy's past 
star formation history (e.g. Mihos & Hernquist 1996; Cox 
et al. 2008). It is also likely to enable supermassive 
black hole growth and the fueling of AGN (e.g.. Heck- 
man et al. 1986; Springel et al. 2005a; Somerville et al. 
2008). Cox et al. (2008) used Smooth Particle Hydro- 
dynamical (SPH) simulations to show that the timescale 
over which merger-induced starbursts arc active depends 
sensitively on the treatment of poorly-understood feed- 
back and ISM physics; they demonstrate that future ob- 
servational constraints on this timescale may provide a 
means to constrain feedback models (Barton et al. 2007, 
and references therein). Historically, SPH simulations 
have treated star forming gas as isothermal, and this 
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treatment results in starburst timescalcs in major merg- 
ers that are quite short-lived, with ~ 100 Myr (e.g. 
Mihos & Hernquist 1996; Cox et al. 2008). Recently, 
it has become popular in SPH simulations to impose a 
stiff equation of state for star-forming gas in order to 
mimic the effects of a multi-phase ISM and to suppress 
star formation and disk fragmentation (Yepes et al. 1997; 
Springel & Hernquist 2003; Governato et al. 2007). Cox 
et al. (2008) showed that a stiff equation of state of this 
kind significantly lengthens the timescale for starburst 
activity in major mergers to ~ 500 Myr. Below we in- 
vestigate the evolution of merger fractions with 100 Myr 
and 500 Myr as a first-order means of addressing the dif- 
ferences between merger- induced starburst fractions in 
different feedback schemes. 

A second obscrvationally-relevant consequence of 
mergers is morphological disturbance. Very large merg- 
ers, especially those with moderately low gas fractions, 
likely play a role in transforming late type disk galaxies 
into ellipticals (e.g., Toomre & Toomre 1972; Barnes & 
Hernquist 1996; Robertson et al. 2006a, b; Burkert et al. 
2008). If gas fractions are high in major mergers (as 
expected at high redshift) then they may play a role 
in building early disks (Robertson et al. 2006a; Hopkins 
et al. 2009, 2008; Robertson & Bullock 2008). More com- 
mon are moderate-size (m/M > 0.1) dark matter halo 
mergers (Stewart et al. 2008), which can produce mor- 
phological signatures like disk flaring, disk thickening, 
and ring and bar-like structures in disk galaxies (Barnes 
& Hernquist 1996; Kazantzidis et al. 2008; Younger et al. 
2007; Villalobos & Helmi 2008; Purcell et al. 2009) as 
well as tidal features seen in massive elliptical galaxies 
(Feldmann et al. 2008). 

Below we explore two possibilities for the evolution of 
the morphological relaxation time with redshift. First, 
we explore a case where the remnant relaxation time 
scales with redshift, approximated by the dark matter 
halo dynamical time (r oc (1 -f z)~", a ~ 1.1 — 1.5; see 
below) , and second we investigate the possibility that re- 
laxation times remain constant with redshift at t ~ 500 
Myr. The latter timescale is motivated by the results of 
Lotz et al. (2008b) who studied outputs from SPH merger 
simulations of z = galaxies in great detail (see Cox et al. 
2006; Jonsson et al. 2006; Rocha et al. 2008; Cox et al. 
2008, for additional descriptions of these simulations and 
their analysis). These choices bracket reasonable expec- 
tations and allow us to provide first-order estimates for 
the evolution in the morphologically disturbed fraction 
with redshift. More simulation work is needed to deter- 
mine how the relaxation times of galaxy mergers should 
evolve with redshift, including an allowance for the evo- 
lution in approach speeds, galaxy densities, and orbital 
parameters (if any). 

Though not discussed in detail here, a third conse- 
quence of mergers is the direct, cumulative deposition 
of cold baryons (gas and stars) into galaxies. For this 
question, one is interested in the full merger history of 
individual objects, rather than the instantaneous merger 
rate or recent merger fraction. Specifically, one may ask 
about the total mass that has been deposited by major 
mergers over a galaxy's history. We focus on this issue 
in Stewart et al. (2009). 

In what follows we use a high-resolution dissipation- 
less cosmological LCDM N-body simulation to investi- 



gate the merger rates and integrated merger fractions of 
galaxy dark matter halos of mass M = 10^^ — 10^^h~^ Mq 
from redshift z = to 4. We adopt the simple technique 
of monotonic abundance-matching in order to associate 
dark matter halos with galaxies of a given luminosity 
or stellar mass (e.g., Kravtsov et al. 2004; Conroy et al. 
2006; Berrier et al. 2006; Conroy & Wechsler 2009), and 
make predictions for the evolution of the galaxy merger 
rate with redshift. 

The outline of this paper is as follows. In §2 we discuss 
the numerical simulation used and the method of merger 
tree construction, while we present merger statistics for 
dark matter halos in §3. In §4 we discuss the method of 
assigning galaxies to dark matter halos both as a function 
of stellar mass, and alternatively as a function of galaxy 
luminosity (compared to L^:(z)). In §5 we present our 
principle results, which characterize the merger rate of 
galaxies as a function of redshift, with comparison to 
observed properties of bright galaxies. We summarize 
our main conclusions in §6. 

2. SIMULATION 

We use a simulation containing 512^ particles, each 
with mass nip = 3.16 x IO^H'^Mq, evolved within a 
comoving cubic volume of 80h~^ Mpc on a side us- 
ing the Adaptive Refinement Tree (ART) A^-body code 
(Kravtsov et al. 1997, 2004). The simulation uses a flat, 
ACDM cosmology with parameters f^M = 1 — = 0.3, 
h = 0.7, and CTg = 0.9. The simulation root computa- 
tional grid consists of 512'^ cells, which are adaptively 
refined to a maximum of eight levels, resulting in a peak 
spatial resolution of 1.2h~^ kpc (comoving). Here we 
give a brief overview of the simulation and methods used 
to construct the merger trees. They have been discussed 
elsewhere in greater detail (AUgood et al. 2006; Wechsler 
et al. 2006; Stewart et al. 2008) and we refer the reader 
to those papers for a more complete discussion. 

Field dark matter halos and subhalos are identified us- 
ing a variant of the bound density maxima algorithm 
(Klypin et al. 1999). A suhhalo is defined as a dark mat- 
ter halo whose center is positioned within the virial ra- 
dius of a more massive halo. Conversely, a field halo is 
a dark matter halo that does not lie within the virial ra- 
dius of a larger halo. The virial radius R and mass M are 
defined such that the average mass density within R is 
equal to Avir (— 337 at z = 0) times the mean density of 
the universe at that redshift. Our halo catalogs are com- 
plete to a minimum halo mass of M ~ IO^^H'^^Mq, and 
our halo sample includes, for example, ^ 15, 000(10, 000) 
and 2,000(500) field halos &t z ~ 0(3) in the mass bins 
10"-^2 and IO^^-^^/i-^Mq, respectively. 

We use the same merger trees described in Stewart 
et al. (2008), constructed using the techniques described 
in Wechsler et al. (2002) and Wechsler et al. (2006). 
Our algorithm uses 48 stored timesteps that are ap- 
proximately equally spaced in expansion factor between 
a = (1 + z)-^ = 1.0 and a = 0.0443. We use standard 
terminologies for progenitor and descendant. Any halo 
at any timestep may have any number of progenitors, 
but a halo may have a single descendant — defined to be 
the halo in the next timestep that contains the majority 
of this halo's mass. The term main progenitor is used 
to reference the most massive progenitor of a given halo, 
tracked back in time. 
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Throughout this work we present results in terms of 
the merger ratio of an infaUing object, m/M ^ where wc 
always define m as the mass of the smaller object just 
prior to the merger and M is the mass main progenitor 
of the larger object at the same epoch. Specifically, M in 
the ratio does not incorporate the mass m and therefore 
m/M has a maximum value of 1.0. Except when ex- 
plicitly stated otherwise, we always use dark matter halo 
masses to define the merger ratio of any given merger 
event, and we always define the merger ratio as the mass 
ratio just before the smaller halo falls into the virial ra- 
dius of the larger one. Because there is not a simple 
linear relation between halo mass and galaxy stellar (or 
baryonic) mass, this is an important distinction. For ex- 
ample, our major mergers, defined by halo mass ratios, 
may not always correspond to major galaxy mergers as 
defined by stellar or baryonic mass ratios (see e.g. Stew- 
art 2009). 

In what follows we investigate two types of mergers. 
The first and most robust of our predicted rates is the 
infall rate: the rate at which infaUing halos become sub- 
halos, as they first fall within the virial radius of the main 
progenitor. These are the results we present in Section 3, 
which describes our 'universal' merger rate function for 
dark matter halos. The second rate is aimed more closely 
at confronting observations and is associated with cen- 
tral mergers between galaxies themselves. Specifically we 
define the destruction rate by counting instances when 
each infaUing subhalo loses 90% of the mass it had prior 
to entering the virial radius of the larger halo. ^ Wc 
are unable to measure central crossings directly because 
the time resolution in our snapshot outputs (typically 
At ~ 250 Myr) is comparable to a galaxy-galaxy cross- 
ing time at the centers of halos, however, for mergers 
with mass ratios > 1/3(1/10), subhalo destruction typ- 
ically takes place ^ 2(3) Gyr after infall to the virial 
radius. Based on simulations of galaxy mergers, this def- 
inition leads to subhalo "destruction" sometime after its 
first pericenter (and likely after second pcricenter), but 
probably before final coalescence (Boylan-Kolchin ct al. 
2008). Note that this second rate (destruction) is more 
uncertain than the first (infall) because, in principle, the 
orbital evolution of infaUing galaxies will depend upon 
the baryonic composition of both the primary and sec- 
ondary objects. Fortunately, as presented in detail below 
(see Table 1) for the relatively high mass-ratio merger 
events we consider, the merger rates (and their evolution 
with redshift) do not depend strongly on whether we de- 
fine a merger to occur at halo infall or at this central 
mass-loss epoch. 

3. DARK MATTER HALO MERGER RATES 

We begin by investigating infall and destruction merger 
rates as a function of mass, merger ratio, and redshift. 
Merger rates are shown for several of these choices in the 
four panels of Figure 1. The upper panels show merger 
rates per unit time for > m/M mass ratio objects falling 
into host halos of mass 10^^ (black lines and crosses) 
and IQ^^h^^MQ (red lines and squares) at three differ- 

^ Subhalo masses are defined to be the mass within a truncation 
radius Rt, which is set to be the minimum of the virial radius and 
the radius where the subhalo density profile begins to encounter 
the background halo density. 



ent redshifts: z ~ Q (solid), z — 2 (dashed), and z ~ i 
(dot-dashed). Host halo mass bins span A log^o M = 0.5, 
centered on the mass value listed. The upper left panel 
presents rates measured at subhalo infall — i.e., the 
merger rates of distinct halos — and the upper right 
panel presents rates of subhalo destruction (when the 
associated subhalo loses 90% of the mass it had prior 
to entering the virial radius of the larger halo), which 
we expect to more closely trace the galaxy merger rates. 
The lower left panel presents infall rates, now plotted 
at a fixed mass ratio (> m/M = 0.1, ...,0.7 from top 
to bottom) and host mass (M = lQ^'^-^h~^ Mq, trian- 
gles; M = lO^^'^/i~"'^M0, squares) as a function of red- 
shift. The same information is presented in the lower 
right panel, but now presented as the rate per unit red- 
shift instead of per unit time. We see that merger rates 
increase with increasing mass and decreasing mass ratio, 
and that the merger rate per unit time increases with 
increasing redshift out to z ^4. 

Wc quantify the measured dependencies using simple 
fitting functions. The merger rate (for both infall and 
destruction rates) per unit time for objects with mass 
ratios larger than m/M into halos of mass M at redshift 
z is fit using 

dN 

^{>m/M) = Mz,M)F{m/M). (1) 

For the infall rate, we find that the normalization 
evolves with halo mass and redshift as At{z,M) = 
0.02 Gyr"^ (1 + z)^ ^ Mfj with M12 the mass in units of 
IQ^^H^^Mq and b = 0.15. The merger mass ratio depen- 
dence is fit by 

(^)'(-S)'. 

with c = 0.5, and d = 1.3. A similar fit describes the 
destroyed rate, as summarized in Table 1. The fits are 
illustrated by solid lines that track the simulation points 
in each of the dN/dt panels in Figure 1. 

The solid and dotted lines in the lower-right panel 
of Figure 1 show that the infall rate per unit redshift, 
dN/dz, is well described by the same mass-dependent 
function, but with a normalization that is only weakly 
dependent on redshift: 

dN 

— {> m/M) = AAz,M)F{m/M), (3) 
az 

where A,{z,M) = 0.27 {dSc/dz f M^^^^ (for infaU rate; 
see Table 1 for destruction rate). As discussed by 
Fakhouri & Ma (2008) (hereafter FM08), a redshift evo- 
lution of this form is motivated by the expectations of 
Extended Prcss-Schcchtcr theory. Note that since dSc/dz 
asymptotes to a constant ~ 1.3 for z > 1 and evolves only 
mildly to ~ 0.9 at z ~ 0, the overall redshift dependence 
is weak. 

To a large extent, our results confirm and agree with 
those of FM08, who studied merger rates for halos in 
the Millennium simulation (Springel et al. 2005b) and 
presented a fitting function for the merger rate per unit 
redshift per unit mass-ratio for halos as a function of 
mass and redshift (the differential of our rate, dN/dz, 
with respect to the merger rate m/M), and concluded 
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Fig. 1. — Dark matter halo infall and destruction rates (see §2) as a function of mass, merger-mass ratio, and redshift. Host halo mass bins 
span Alogj^Q M = 0.5. Top Left: Infall rate per Gyr as a function of merger mass ratio. The dashed (pink) lines are a comparison to the 
results of Fakhouri & Ma (2008) for M = 10^^/i"^Mq (lower) and M = 10^^h~'^MQ (upper) halos. Top Right: Identical to top left, but 
for the destruction rate of halos instead of the infall rate. Bottom Left: Infall rate per Gyr as a function of redshift. Bottom Right: Infall 
rate per unit redshift, as a function of redshift. In the top panels, black and red lines correspond to different host halo masses (10^^ and 
10^^h~^ Mq), while in the bottom panels, magenta, black, blue and red lines correspond to different merger ratios (m/M > O.l, 0.3, 0.5, 0.7). 
In all panels, the mass ratio, m/M is defined prior to the infall of the smaller system. Error bars are Poissonian based on the number of 
host halos and the number of mergers. Horizontal error bars on the bottom figures have been omitted for clarity, but are identical to those 
in Figure 2. 



that it was nearly universal in form. For comparison, the 
pink dashed lines in the top left panel of Figure 1 show 
the implied expectations based on the FM08 fit for M = 
lO^^h-^MQ (lower lines for each pair) M = 10^^h~'^MQ 
(upper lines) halos. ^ The agreement is quite remarkable, 
especially in light of the fact that the simulation, merger 
tree algorithm, and halo finder all differed substantially 
from our own. Note that the agreement is particularly 
good over the mass ratios m/M ^ 0.05 — 0.5, that are 
likely the most important for galaxy formation (in terms 
of their potential for morphological transformation and 

* We use their fit for the 'stitching' merger rate, which corre- 
sponds most closely to our own definition for halo mergers. 



overall mass deposition, see Stewart et al. 2008). We note 
however that our infall rate data are smaller than FM08 
by a factor of ~ 1.5 for very large mass-ratio mergers 
m/M > 0.7 and by a factor of ~ 2 for very small mass- 
ratio mergers m/M < 0.01 (this discrepancy for small 
mergers is slightly worse at low redshift, z < 0.3). In 
addition, we find a slightly stronger mass dependence, 
dN/dt oc as opposed to dN/dt oc M°-i as found 

by FM08. 

It is interesting to note that in an independent analysis 
of the Millennium simulation, Gencl et al. (2009) studied 
the (infall) merger rates of halos by defining halo masses 
and mergers in slightly different ways from FM08, in an 
effort to further remove artifacts of the halo-finding al- 
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TABLE 1 

Merger Rate Fitting Function Parameters for Equations 3.1,3.2, and 3.3. 



Dark Matter Halos: 

(M = lO"-"~i^-^/i~Uf0) 


A(z, A'/)t 


c 


d 


dN/dt: (INFALL) simple fit 


0.020 (1 + zf-'-^ A/O;^!^ 


0.50 


1.30 


dN/dt: (INFALL) complex fit 


0.020 (1 + 2)2-3 A//,y'^ 


0.4+ .05z 


1.30 


dN/dz: (INFALL) 


0.27 {dSc/dz)'^ A/o/s 


0.50 


1.30 


dN/dt: (DESTROYED) 


0.022 (1 + z)'^-'^ Mj";^^ 


0.54 


0.72 


dN/dz: (DESTROYED) 


0.32 (dS^/dz) A^iY 


0.54 


0.72 


Galaxy Luminosity Cuts: 

(L > /L*, O.K / < 1.0) 




c 


d 


dN/dt: 


0.02(1 + /) (1 + 2)^-^ 


0.54 


0.72 


Merger Fraction in past T Gyr: 
(Frac < 0.6, T < 4) 


0.02T(l + /)(l + z)^-0 


0.54 


0.72 


Galaxy Stellar Mass Ranges: 

(F(x) = F(m,/KU)) 


A(2)t 


c.t 


dj 


dN/dt (lOi"-"*/© < h'h < 10"^-^A/o): 


0.015 e^ o^ 


0.30 


1.1 - 0.2z 


dN/dt (1O"'-^M0 < A/* < lO^^-"Af0): 


0.035 e" '''^ 


0.25 


1.1 - 0.22 


dN/dt (lO"-"Af0 < A/*): 


0.070 e^ o^ 


0.20 


1.0 - 0.32 



^When not dimensionless, units are Gyr ^. 

tMass-ratio variable for galaxy stellar mass merger rates are identified with a stellar mass ratio, r = m»/Af*. 



TABLE 2 

Dark Matter Halo Mass-Luminosity Relationship by Number Density Matching. 



z 


Source & Rest-Frame Band 


> O.IL, 


> 0.4L, 


> L, 


r{z) 
[Gyr] 


Ug T 


Mdm 


Ug T 


Mdm ^ 


Ug T 


Mdm ^ 


0.1 


SDSS (rO-i-band) ^ 


29 


1011.2 


10 


lOil.V 


3.2 


1012.3 


1.79 


0.3 


COMBO-17/DEEP2 (B-band) ^ 


20 


1011.4 


5.8 


1012.0 


1.5 


1012.6 


1.50 


0.5 


COMBO-17/DEEP2 (B-band) ^ 


24 


1011.3 


7.0 


1011.9 


1.8 


1012.5 


1.28 


0.7 


COMBO-17/DEEP2 (B-band) ^ 


20 


1011.4 


5.8 


1012.0 


1.5 


1012.6 


1.09 


0.9 


COMBO-17/DEEP2 (B-band) ^ 


25 


1011.3 


7.3 


1011.9 


1.9 


10l2.,S 


0.95 


1.1 


COMBO-17/DEEP2 (B-band) ^ 


19 


1011.4 


5.5 


1012.0 


1.4 


1012 = 


0.81 


2.2 


Keck Deep Fields (UV) 


20 


1011.3 


6.4 


1011.7 


1.8 


1012.2 


0.49 


~ 2.6 


Extrapolation (UV) 


~ 18 


1011.3 


~ 5.2 


lOii.s 


~ 1.2 


1012.3 


0.40 


3 


Keck Deep Fields (UV) ^ 


15 


1011.2 


3.8 


1011.7 


0.90 


1012.2 


0.32 


3 


Lyman Break Galaxies (V-band) * 


NA 


NA 


5.0 


1011.7 


0.82 


1012.2 


0.32 


4 


HUDF + HST ACS Fields (UV) ^ 


18 


10"" 


3.2 


1011.6 


0.61 


1012.0 


0.25 



tlO-3 Mpc-3 
iBlanton et al. (2003) 

^Faber et al. (2007) with DEEP2 optimal weights. 
^Sawicki & Thompson (2006) 

*Shapley et al. (2001) - note that rest-frame V number densities 
match well with > 0.4L* and > L, values in rest-UV at 2 = 3. 
^Bouwens et al. (2007) 

gorithm of the simulation. Among other results, their 
findings suggested that the merger rates from FM08 are 
slightly too high (by < 50%) for low redshift and for mi- 
nor (< 1/10) mergers. This is qualitatively similar to 
the differences between FM08 and our own results, moti- 
vating the need for future study regarding the sensitivity 
of merger statistics from dark matter simulation on halo 
finding algorithms, as well as halo mass and merger def- 
initions. 

Our results also largely agree with an investigation of 
the major merger rate (> 1/3 mergers) of halos and sub- 
halos by Wetzel et al. (2009). They found that the infaU 
rate for halos (M = 10" - lO^^/i-^M©) evolves with 
redshift as dN/dt = ^(1 + z)" (with A - 0.03 and 
a = 2.0 — 2.3) from z = 0.6 — 5, in good agreement with 
our infall rates both in slope and in normalization (see 



Table 1). Wetzel et al. (2009) also reported on the sub- 
halo merger rate in their simulations (the rate at which 
satellite subhalos finally merge with the central subhalo) 
and found similar behavior as field halos for low redshift 
{A ~ .02 and a ~ 2.3 for z = 0.6 — 1.6, in good agree- 
ment with our destruction rates) but with a significantly 
flatter slope for high redshift {A ^ 0.08 and a = 1.1 for 
z = 2.5 — 5, a factor of 2 — 3 lower than our results, with a 
significantly flatter slope). Even though our destruction 
rate attempts to track a similar physical phenomenon as 
their subhalo merger rate — the rate of impact of satellite 
galaxies onto central galaxies — we find destruction rates 
to show qualitatively similar behavior to infall rates at 
all redshifts. 

Wc speculate that the discrepancy between their re- 
sults and ours may be due primarily to differences in 
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definition. For example, we define an infalling halo to be 
"destroyed" once it loses 90% of its infall mass (see §2). 
Wetzel et al. (2009) defines subhalo mergers by track- 
ing the evolution of the subhalo's 20 most-bound parti- 
cles, resulting in a much more stringent definition of a 
merger, and increasing the time delay between infall to 
the virial radius (iinf ) and the time at which the satellite 
is destroyed (imcrgc)- More importantly, we define the 
merger mass ratio by the halo masses when the satellite 
halo first falls into the virial radius of the host, mi„f/Minf . 
Although we track the subhalo until it has lost 90% of 
its mass in order to assign a proper time that the merger 
takes place, we do not redefine this merger ratio based on 
any subsequent growth or decay of either halo. Although 
Wetzel et al. (2009) defines the satellite halo's mass in an 
identical fashion, they allow for the growth of the cen- 
tral halo during the decay time of the subhalo. Once the 
subhalo is destroyed, they use the host halo mass at this 
time (minus the mass of the subhalo, so that m/M < 1) 
and thus define the merger ratio as minf/Mmcrgc- As a 
consequence, the host halo has a significant time period 
(Emerge" ^inf) to grow in mass, leading to smaller mass ra- 
tio definitions for identical merger events, as compared to 
our definition. This effect is likely negligible at late times, 
when halos do not grow significantly over the ~ 2 Gyr 
decay timescales typical for major mergers. This may be 
why the two studies agree rather well for z < 1.6. How- 
ever, the central halo's mass growth on these timescales 
becomes increasingly important at high redshift, possibly 
explaining the flattening of a reported by Wetzel et al. 
(2009), as compared to our own results. 

Fakhouri & Ma (2009) investigates this issue to some 
degree by studying the subhalo merger rate in the Mil- 
lennium simulation using differing mass ratio definitions. 
Whether they implement a merger ratio definition simi- 
lar to our destruction rate, or one more similar to that 
of Wetzel et al. (2009), their merger rates remain well-fit 
to a power law in (1 -I- z), in line with our results. In this 
case, the underlying cause of the discrepancy between 
our merger rates (and those of Fakhouri & Ma (2009)) 
and the subhalo merger rates reported by Wetzel et al. 
(2009) remains uncertain. Such comparisons between 
our respective results highlight the differences that are 
manifest in defining mergers. When including baryons, 
properly defining mergers and merger mass ratios be- 
comes even more complicated (see e.g. Stewart 2009), but 
even between dark matter structures, differences such as 
those found between our work, Wetzel et al. (2009) and 
Fakhouri & Ma (2009) further motivate the need for fo- 
cused simulations in order to determine the timescales 
and observational consequences associated with the much 
more cleanly defined rate with which dark matter sub- 
halos first fall within the virial radii of their hosts. 

4. ASSOCIATING HALOS WITH GALAXIES 

While dark matter halo merger rates at a given mass 
are theoretically robust quantities to compute in our 
simulation, they are difficult to compare directly with 
observations. One particularly simple, yet surprisingly 
successful approach is to assume a monotonic mapping 
between dark matter halo mass M (or similarly the 
halo maximum circular velocity) and galaxy luminosity 
L (Kravtsov et al. 2004; Tasitsiomi et al. 2004; Vale & 
Ostriker 2004; Conroy et al. 2006; Berrier et al. 2006; 



PurceU et al. 2007; Marfn et al. 2008; Conroy & Wech- 
sler 2009). With this assumption, provided that we know 
the cumulative number density of galaxies brighter than 
a given luminosity, ng(> L), we may determine the asso- 
ciated halo population by finding the mass M{L) above 
which the number density of halos (including subhalos) 
matches that of the galaxy population nh{> Mdm) = 
ng{> L). Table 2 shows the number densities of vari- 
ous galaxy populations from redshifts z = 0.1 — 4 ob- 
tained using a variety of surveys for galaxies brighter 
than L ^ fL^, where / = 0.1,0.4, and 1.0. We list 
the associated number-density matched minimum dark 
matter halo mass in each case, Mdm, and we use this 
association to identify halos with galaxies below. For ex- 
ample, from the top left entry of this table, we see that 
nh{> 10"-2/i-iMo) = ng{> O.IL*) at z = 0.1. 

One important point of caution is that the luminos- 
ity functions used to make these assignments at different 
redshifts vary in rest-frame band, as indicated in Column 
2. Specifically, one concern might be that UV luminosity 
at low redshift is not strongly correlated with dark mat- 
ter halo mass, so assuming such a correlation for high 
redshift galaxies is not valid. Unlike their low redshift 
counterparts, however, there is a strong correlation in 
high redshift (z > 2) galaxies between star formation and 
total baryonic mass, as well as a trend for more UV lumi- 
nous galaxies to be more strongly clustered, suggesting 
that connecting UV luminosity to halo mass at these red- 
shifts is a valid technique (see discussion in Conroy et al. 
2008, and references therein). Encouragingly, as shown 
in the two z = 3 rows, the number density of / galax- 
ies from Sawicki & Thompson (2006; rest-frame UV) and 
Shapley et al. (2001; rest- frame V) are quite similar. 

We also note that the data from these various sources 
will contain uncertainties in the number counts of 
galaxies from, e.g. cosmic variance. For example, the 
COMB017/DEEP2 data fluctuates about a nearly con- 
stant value (- 0.02 - O.OS/i^Mpc-^) from z = 0.3 - 1.1, 
suggesting a ~ 30% uncertainty in these values. We 
find that a 30% error in the observed number density 
typically translates into a similar 30% error in the as- 
signed minimum halo mass in our simulation. Since 
dark matter halo merger rates are only weakly depen- 
dent on halo mass (oc M^^), this should result in only a 
~ 10% uncertainty in our merger rates. Thus, the merger 
rates we present here should be relatively robust to small 
errors in observational uncertainties. For example, if 
we adopt minimum halo masses (regardless of redshift) 
of Mdm = 10"■^ lO"-'^, IO^^-^H-^Mq as corresponding 
to > 0.1, 0.4, l.Oi* galaxies, respectively, our resulting 
merger rates change by < 25% (typically 5 — 15%). 

A related approach is to use observationally-derived 
stellar mass functions and to assume a monotonic rela- 
tionship between halo mass and stellar mass M*. Though 
a monotonic relationship between total stellar mass and 
dark matter mass avoids the issue of color band that 
arises in luminosity mapping, we cannot use it explore 
merger rates as a function of stellar mass above z ~ 2 
because the stellar mass function is poorly constrained 
beyond moderate redshifts. For our analysis, we will 
adopt the relation advocated by Conroy & Wechsler 
(2009, hereafter CW09; interpolated from the data shown 
in their Figure 2). For example, CW09 find that the 
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Fig. 2. — Expected merger rates per Gyr for galaxies of an indicated type as a function of redshift. The vortical error bars show Poisson 
errors on both the number of main halos and the total number of mergers averaged over per redshift bin while the horizontal error bars 
show the redshift bins used to compute the merger rate at each redshift. The error bars do not include uncertainties in the mapping of 
mass to luminosity or stellar mass. Left: Merger rate into galaxies with L > Lt, involving objects with total mass ratios m/M > 0.1, ...,0.7 
as indicated. Middle: Merger rate into galaxies with L > 0.1L«. Right: Merger rates for galaxies of a given stellar mass involving objects 
with stellar mass ratios mt/Mt, > 0.3 and 0.6 as a function of redshift. (Note that the redshift range in this panel only goes to 2: = 2.) 
We include two different stellar mass cuts in this panel, represented by the red and blue lines. The dotted lines in this panel show an 
extrapolation out to z ~ 4, based on our fit to the z < 2 simulation data. The filled diamonds show observational results for the same 
stellar mass cuts from Bundy et al. (2009). 



halo mass M associated with stellar masses of 7\jf, = 
(1,3,10) X IQi^Mg at z = 0,1,2 are M{z = 0) ~ 
(2.5,7.0,47) X 10"/i-^Mo; M(z = 1) = (4.0,9.6,41) x 
10"; and M{z = 2) = (2.1,3.9,10) x lO^^h-^MQ. We 
note that because this mapping between stellar mass 
and halo mass is not well fit by a constant ratio, Af* = 
/Mdm, merger rates in terms of stellar mass ratios show 
qualitatively different evolution with redshift (see §5.1). 
This is primarily because mergers of a fixed dark mat- 
ter mass ratio do not typically correspond to the same 
stellar mass ratio (see Stewart 2009). 

Note that while the dark matter halo merger rates 
presented in §3 give robust theoretical predictions, the 
merger rates we will present in terms of luminosity (or 
stellar mass) are sensitive to these mappings between 
halo mass and L (or M*). In addition, it is difficult to 
perform a detailed investigation into the errors associated 
with these mappings, as there arc inherent uncertainties 
in the luminosity and stellar mass functions, especially 
at z > I. It is also possible that the monotonic map- 
ping between halo mass and L (or M*) may break down 
at z > 1 (see discussion in CW09). These uncertain- 
ties must be kept in mind when comparing our predicted 
merger rates (in terms of L or Af*) to observations, espe- 
cially at high redshift. Nevertheless, the halo masses we 
have associated with a given relative brightness should 
be indicative. 

5. GALAXY MERGER PREDICTIONS 

5.1. Merger Rates 

Our predicted merger rates (per galaxy, per Gyr) and 
their evolution with redshift, averaged over L > L^, 
and L > 0.1 galaxy populations, are illustrated in 
the left and middle panels of Figure 2. Rates are pre- 
sented for a few selected dark matter halo mass ratio 
cuts m/M > 0.1, 0.3, 0.5, and 0.7. Here, galaxy merger 
rates are defined using the destruction rate, when the 
infalling subhalo is identified as destroyed in the simula- 
tion (see §2). The solid lines correspond to a fit in the 
form of Equation (1), with the normalization evolving 
as At{z,f) (X (1 -I- /) (1 -I- z)"^-^ for L > fL^, galaxies. 



The explicit best-fit parameters for the merger rate as a 
function of luminosity cut are given in Table 1 . 

For comparison, the right panel in Figure 2 shows the 
predicted evolution in the merger rates per galaxy for 
two bins of stellar mass, according to the CW09 map- 
ping described above: lO^^^Afo < A/, < 10^°-^Mq 
(lower, blue) and A/, > IO^^Mq (upper, red). Shown are 
merger rates for two choices of stellar mass ratio merg- 
ers, (m*/Af*) > 0.3, 0.6 (solid and dashed lines respec- 
tively). The solid and dashed lines correspond to fits to 
our simulation results in the form of Equation (1), with 
At{z) oc e^'^^. The explicit best-fit parameters for these 
two stellar mass bins (as well as an intermediate bin, 
10^°-^Mq < at* < lO" °A/0) can be found in Tabic 1. 
Table 1 also provides best fit parameters for the function 
F{r) (Equation 2) where now we associate the ratio r 
with the stellar mass ratio r = m^,/M^,. Note that we 
only show our simulation points for z < 2 in this panel, 
due to uncertainties in the stellar mass function at high 
redshift. 

While the galaxy merger rate cannot be observed di- 
rectly, it can be inferred using a number of different 
techniques. Mergers that are about to occur may be 
forecast by counting galaxy close pairs, and close pair 
fractions are often used as a proxy for the merger rate. 
The filled diamonds in the right panel of Figure 2 are re- 
cent merger-rate estimates from the pair count study of 
Bundy et al. (2009), for the same two stellar mass bins 
shown in the simulations (blue for the lower mass bin, 
red for the upper mass bin). Bundy et al. (2009) have 
used the simulation results of Kitzbichler & White (2008) 
to derive merger rates from the observed pair fraction. 
Overall, the trends with mass and redshift are quite sim- 
ilar and this is encouraging. However, the Bundy et al. 
(2009) results correspond to mergers with stellar mass- 
ratios larger than m*/A/* > 0.25. Our normalization 
is a factor of ~ 2 too high compared to this, and only 
matches if we use larger merger-ratios m^/M^ > 0.5. It 
is possible that this mismatch is associated with the dif- 
ficulty in assigning merger timcscalcs to projected pairs 
(see, e.g. Berrier et al. 2006). It may also be traced 
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back to uncertainties in assigning stellar masses to dark 
matter halo masses, however, since merger rates have 
relatively weak dependence on halo mass, it would re- 
quire increasing our assigned stellar masses by a factor 
of 3 in order to account for this discrepancy solely 
by errors in assigning stellar mass (such an increase in 
stellar mass would result in unphysical baryonic content 
for dark matter halos: e.g., 10^^/i~^Mq halo containing 
Ah > lO^/i-^Mo). 

There are a number of other observational estimates 
of the merger rate based on pair counts of galaxies (e.g., 
Patton et al. 2002; Lin et al. 2004; Bell et al. 2006; Kar- 
taltcpe et al. 2007; Kampczyk et al. 2007; de Ravel et al. 
2009; Lin et al. 2008; Mcintosh et al. 2008; Patton & At- 
field 2008; Ryan et al. 2008). We choose to compare our 
results to Bundy ct al. (2009) as a recent representative 
of such work; primarily because it is more straightfor- 
ward for us to compare to samples that are defined at 
a fixed stellar mass and stellar mass ratio. It is also 
difficult to compare to many different observational re- 
sults on the same figure self-consistently, because differ- 
ent groups adopt slightly different cuts on stellar mass (or 
luminosity) and on mass ratios for pairs. We note that if 
we were to extrapolate our best-fit curves to higher red- 
shift than our data (2; ~ 4), we find good agreement be- 
tween our simulation data and the merger rate estimates 
using CAS (concentration, asymmetry, dumpiness) mor- 
phological classifications from Consclicc ct al. (2003) for 
galaxies with > 10^^h~^ Mq. However, the mapping 
between stellar mass and halo mass adopted from Con- 
roy & Wechsler (2009) is only valid to 2; = 2 (and most 
robust for z < 1), so extrapolating these fits to 2; ~ 4 is 
only a first-order check, and should not be considered a 
reliable prediction. 

5.2. Merger Fractions 

Another approach in measuring galaxy merger rates 
is to count galaxies that show observational signatures 
of past merging events such as enhanced star formation, 
AGN activity, and morphological disturbances. Unfortu- 
nately, the timcscale over which any individual signature 
will be observable is often extremely uncertain, and will 
depend on the total mass and baryonic makeup of the 
galaxies involved as well as many uncertain aspects of 
the physics of galaxy formation (e.g. Berrier et al. 2006; 
Cox et al. 2008; Lotz et al. 2008b; Kitzbichler & White 
2008). In order to avoid these uncertainties, we present 
results for merger fractions using several choices for look- 
back timescale here. 

The three panels of Figure 3 show the predicted evo- 
lution of the merger fraction in galaxies brighter than 
0.4L* for three different choices of merger lookback time 
and for various choices for the total mass merger frac- 
tion m/M > 0.1,. ..,0.7. The horizontal error bars on 
this figure show the actual rcdshift bins used to compute 
the merger fractions. The left and right panels show 
the merger fraction within 100 Myr and 500 Myr, re- 
spectively ^, and the middle panel shows the fraction of 
galaxies that have had a merger within the past halo dy- 

^ In most cases, the available timesteps (At ~ 250 Myr) arc too 
widely spaced to directly measure fractions within 100 Myr. For 
this reason, the left panel is actually the merger fraction within the 
last At timestep, scaled down by a factor of (At/lOOMyr) ~ 2.5. 



namical time ^ t{z), where t(z) ~ 2.0 Gyr (1 -I- z) ^-^^ 
for z<l and r(z) ~ 2.6 Gyr (1 -I- z)''^-^ for z > 1. 

5.2.1. Merger-driven starbursts 

Several recent studies of star formation rates in galax- 
ies at 2; = — 1 suggest that the cosmic SFR density is 
not dominated by strongly disturbed systems with brief 
periods of intense star formation, as might be expected 
if merger-driven starbursts arc common. Instead, the 
SFR density appears to be dominated by normal, non- 
merging galaxies (Wolf et al. 2005; Bell et al. 2005; Jo- 
gee et. al 2008; Noeske et al. 2007). That is, < 30% 
of the instantaneous SFR density at a given rcdshift 
(from z = — 1) is derived from morphologically dis- 
turbed galaxies, which may be currently undergoing a 
merger-induced starburst. Even at high redshift (z ~ 2), 
a comparison of the clustering of star-forming galaxies 
to that of dark matter halos suggests that these galaxies 
are consistent with massive galaxies (in massive DM ha- 
los) quiescently forming stars, as opposed to less massive 
galaxies (less massive DM halos) in the midst of merger- 
induced starbursts (Conroy et al. 2008). However, this 
conclusion is based on the assumption that UV-bright 
galaxies at this redshift comprise a representative sam- 
ple of star-forming galaxies. 

As discussed in the introduction, the briefest timcscalcs 
we expect for merger-triggered starbursts is 100 Myr 
(Mihos & Hernquist 1996; Cox et al. 2008), and for these 
models we expect the SFR to increase to ~ 20 times the 
isolated value for m/M > 0.3 events (Cox et al. 2008). 
(While we adopt these timescales as "typical" of galaxy 
mergers, it is important to keep in mind that Cox et al. 
(2008) focuses on z = galaxies. High redshift galax- 
ies should typically contain higher gas fractions, which 
may impact the properties of merger-induced starbursts 
at these epochs.) As we see from the left-panel of Fig- 
ure 3, the fraction of galaxies that have a merger large 
enough (m/M > 0.3) to trigger such a burst is quite 
small, < 1% for z < 1. It is therefore not surprising 
that stochastic starbursts of this kind do not dominate 
the SFR density at moderate to low redshifts. Even at 
at higher redshift (z = 3 — 4), the fraction of galax- 
ies with major mergers on these timescales is less than 
~ 6% of the total bright galaxy population (consistent 
with the results presented in Somcrville ct al. (2008), for 
their semi-analytic model). However, galaxy gas frac- 
tions are expected to increase with redshift (Erb et al. 
2006), which could presumably result in significant star- 
burst activity from more minor mergers (as well as pro- 
viding fresh gas accretion in a more cumulative sense, 
see Stewart et al. 2009). A higher fraction of galaxies 
have experienced such minor (> 1/10) mergers on these 
timescales at z = 3 — 4 (~ 15%). 

Alternatively, if merger-driven starbursts remain ac- 
tive for 500 Myr, as other models suggest, then their 
enhancements are expected to be less pronounced (with 
an SFR ^ 5 times isolated; Cox et al. 2008). In this 
case, the right panel of Figure 3 is the relevant predic- 
tion, and we see that (at most) ~ 3 — 9% of bright galax- 
ies could exhibit signs of such elevated SFR activity be- 
tween z = and z = 1. It seems that in either case, we 

S We use r = R/V oc (Av(z) puiz))-'^^^ , such that the halo 
dynamical time is independent of halo mass. 



Fig. 3. — The fraction of halos that experience at least one merger larger than m/M in the past 100 Myr (left), halo dynamical time r 
(middle), or 500 Myr (right), as a function of z. Error bars show the Poisson y/N error based on the both the number of main halos and 
the total number of mergers averaged over, while the horizontal error bars show the redshift bins used to compute the merger rate at each 
redshift. The error bars do not include uncertainties in the mapping of mass to luminosity. The symbols represent estimates of the observed 
merger fraction at various rcdshifts, based on Jogee et. al (2008, red crosses) and Lotz ct al. (2008a, green plus signs), respectively. 



would not expect merger-triggered activity to play a ma- 
jor role in driving the integrated star formation rate at 
these epochs. Only at the highest rcdshifts z ^ 3 would 
this seem possible. However, we once again point out 
that the detailed study of Cox et al. 2008, which we have 
quoted here, focuses on low redshift galaxies, with gas 
fractions < 30%. If minor mergers with very high gas 
fractions (> 50%) are capable of triggering starbursts, 
then over half of all bright galaxies at z > 2 (where such 
high gas fractions are more common) may be in the pro- 
cess of starbursting. 

Under the presumption that only major mergers trig- 
ger starbursts, we note that our numbers are an upper 
limit on the fraction of bright galaxies that could be expe- 
riencing merger-induced starbursts, because moderately 
high gas fractions arc also necessary. For example, a 
study of 216 galaxies at z 2 — 3 by Law ct al. (2007b) 
found that galaxy morphology (in rest-frame UV) was 
not necessarily correlated with star formation rate, and 
in a recent examination of two Chandra Deep Field South 
sources using adaptive optics, Melbourne et al. (2005) 
found an example of a merger of two evolved stellar pop- 
ulations, in which the major merger signature was not 
accompanied by a burst of star formation, presumably 
because both galaxies were gas-poor. Lin et al. (2008) 
finds that 8% (25%) of mergers at z - 1.1 (0.1) are 
gas-poor, suggesting that this issue, while less dominant 
at high redshift. is a significant effect and must be taken 
into consideration. 

We note that we have focused on galaxies which are 
in the midst of a merger- induced starburst. The linger- 
ing impact these bursts will have on the cumulative star 
formation histories (SFH) of galaxies in a separate is- 
sue entirely. A recent study by Cowie & Barger (2008) 
traced recent star formation in > 2000 galaxies from 
z = .05 — 1.5 and found that roughly a quarter of these 
galaxies showed color indications (AB3400-AB8140 vs. 
EW in H/3) indicative of starbursts in the past 0.3 — 1.0 
Gyr. Once again, we find our predictions to be broadly 
consistent with this result, with ^ 20% of bright galax- 
ies (> 0.4i*) having experienced a major merger in the 
past Gyr at 2; ~ 1. With this study of individual galax- 
ies' star formation histories emphasizing the importance 
of starbursts, and the previously mentioned studies of 
the global SFR density emphasizing the importance of 



star formation in normal, non-merging systems, we find 
that our predicted merger rates are broadly consistent 
with both results, suggesting that while starbursts may 
not be the globally dominant form of star formation in 
the Universe, they still play an important role in the star 
formation histories of galaxies. Detailed progress in un- 
derstanding the full importance of merger-induced star- 
bursts on the global SFR density of the Universe will re- 
quire a better understanding of the timescales and signa- 
tures associated with galaxy mergers and merger-induced 
starbursts. 

5.2.2. Morphological signatures 

Even if the contribution to the overall SFR due to very 
recent mergers remains low, this does not necessarily im- 
ply that there would be a lack of morphological signa- 
ture. The timescale for morphological relaxation may be 
significantly longer than starburst activity. Though the 
precise timescales for relaxation are uncertain, the mid- 
dle and right panels of Figure 3 explore merger fractions 
for two reasonable choices: a fixed 500 Myr timescale and 
a redshift-dcpcndent halo dynamical time r. 

Lotz ct al. (2008a) used AEGIS survey data to study 
the morphological evolution and implied galaxy merger 
fraction from redshift z = 0.2 to 1.2. The merger fraction 
results for > 0.4L* galaxies from Lotz et al. (2008a) are 
shown by the green pluses in the middle and right pan- 
els of Figure 3. In a similar investigation, Jogee et. al 
(2008) study z ~ 0.2 — 0.8 galaxies using a combination 
of HST, ACS, Combo-17, and Spitzer 24 /im data to esti- 
mate the fraction of "strongly disturbed" galaxies. Their 
results are shown by the red crosses in Figure 3, and 
are in reasonably good agreement with the points from 
Lotz et al. (2008a). We note that the data from these 
very recent works seem to be in good agreement with 
the m/M > 0.3 merger fraction if the relaxation time is 
close to T. This case in particular has a fairly weak evo- 
lution because r is decreasing with time. Interestingly, 
however, due to the rather large measurement uncertain- 
ties, the data are also in reasonable agreement with the 

^ The data from Jogee et. al (2008) correspond to a fixed stellar 
mass cut at M« ~ 2.5 X 10^'' Mq, but the associated dark matter 
halo mappings from CW09 are close to those for galaxies with 
> 0.4Lt (see Table 1 and our discussion in §4). 
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Fig. 4. — The fraction > 0.4L» galaxies at 2 ~ 0, 1, 2, 3 (black solid to red dashed lines) that have experienced a major merger (m/M > 0.3) 
over a given time period. Symbols show the simulation data, while the lines are given by the fit in Table 1. Left: Merger fraction since t, 
normalized by the halo dynamical time at each redshift, t{z = 0,1,2,3) ~ 1.95,0.92,0.49,0.32 Gyr. Right: Merger fraction in the past t 
Gyr. Error bars show the Poisson \/]V error based on the total number of mergers, and are comparable to the symbol sizes. The error bars 
do not take uncertainties in the mapping of mass to luminosity into account. 



fixed relaxation timescale case of 500 Myr (which has a 
steeper evolution with z), as long as more minor merg- 
ers {m/M > 0.1) can trigger the observed activity. The 
fact that the data matches both the predictions in the 
middle panel and right panel of Figure 3 draws attention 
to the inherent degeneracies in this comparison. The 
same merger fractions are obtained with high-mass ra- 
tio merger events and lookback times or with lower mass 
ratio mergers with slightly shorter lookback times. 

We may also compare our predictions with the results 
of Melbourne et al. (2008), who imaged lb z ^ 0.8 lu- 
minous infrared galaxies (LIRGs) with the Keck Laser 
Guide Star (LGS) AO facility, and found that 3/15 of 
the galaxies showed evidence for a minor merger, while 
only 1/15 was consistent with a major merger. These re- 
sults match our expectations for major {m/M > 0.3) and 
minor {m/M > 0.1) merger fractions at z 0.8 fairly 
well, considering the small number statistics. Similarly, 
Shapiro et al. (2008) study 11 rest-frame UV/optical- 
selected z^2 galaxies with spectroscopic data from SIN- 
FONI on the VLT, and estimate that '-^25% of these sys- 
tems are likely undergoing a major (mass ratio < 3:1) 
merger. Again, our expectations as shown in the middle 
and right panels of Figure 3 are consistent with these 
numbers. 

The above discussion makes it clear that meaning- 
ful comparisons between observed morphologically dis- 
turbed fractions and predicted merger fractions rely fun- 
damentally on understanding how the mass ratio in- 
volved affects the morphological indicator and on the 
associated relaxation timescales of the associated rem- 
nants. In addition, merger rates are expected to depend 
sensitively on the galaxy luminosity and redshift (see Ta- 
ble 1). Comparisons between observational results and 
theory therefore require great care, especially as it con- 
cerns the evolution of the merger rate. If, for exam- 
ple, higher redshift measurements are biased to contain 
brighter galaxies than lower redshift measurements, then 



the redshift evolution will likely be steeper than the un- 
derling halo merger rate at fixed mass. Or, if higher red- 
shift measurements are sensitive to only the most massive 
mergers, while lower redshift measurements detect more 
subtle effects, then the evolution in the merger rate will 
be biased accordingly. 

5.2.3. High Redshift Expectations 

As seen clearly in Figures 2 and 3, the merger rate 
per galaxy and the corresponding merger fraction at a 
fixed time are expected to rise steadily towards high red- 
shift. Even after normalizing by the halo dynamical time, 
which decreases with redshift, this evolution with red- 
shift persists, as seen in Figure 3 (middle). This point 
is emphasized in Figure 4, which shows the fraction of 
L > OALf galaxies that have had a merger larger than 
m/M = 0.3 within the last t Gyr (right) and within the 
last t/T{z) (left). The left-hand panel scales out the evo- 
lution in the halo dynamical time. We see that ~ 50% of 
z=3 galaxies arc expected to have had a major merger in 
the last 700 Myr, and that these galaxies are ~ 4 times 
as likely to have had a significant merger in the last dy- 
namical time than bright galaxies at z=0. It would be 
surprising then if mergers did not play an important role 
in setting the the properties of most z — 3 galaxies like 
Lyman Break Galaxies (LBGs). These major mergers 
should (at least) deliver a significant amount of gas to 
fuel star formation, affect LEG dynamics, and perhaps 
trigger starburst activity. If LBGs represent a biased 
sample at z = 3 (of unusually bright galaxies, more likely 
to have recently undergone a merger-induced starburst) ^ 
then it may be possible that the merger fraction in LBGs 
is even higher than the global merger fraction for > 0.4L* 
galaxies. 

At higher rcdshifts, z > 3, wc expect major mergers 

* It is estimated that ~ 75% of all bright galaxies at 2 3 are 
LBGs (Marchesini et al. 2007; Quadri et al. 2007). 
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to become increasingly common. The brightest galax- 
ies L > 0.4L, should be undergoing mergers frequently, 
with an overwhelming majority of z = 4 galaxies having 
experienced some significant merger activity in the last 
~ 500 Myr. 

6. CONCLUSION 

We have used a high-resolution ACDM iV-body sim- 
ulation to investigate the instantaneous merger rate of 
dark matter halos as a function of redshift (from z ~ 
— 4), merger mass ratio, and host halo mass from 
M = 10^^ to l[y^^h~^MQ. Merging companions as small 
as m = IQ^^h'^^MQ were tracked. We use number den- 
sity matching to associate galaxies with dark matter 
halos and present predictions for the merger rate and 
merger fraction as a function of galaxy luminosity and 
stellar mass. The principle goal has been to present raw 
merger statistics that can be compared directly to obser- 
vations of galaxies to high redshift. Fitting functions that 
describe our results as a function of luminosity, mass, 
mass-ratio, and redshift are provided in Table 1. 

Our main results may be summarized as follows: 

1. A simple fitting function describes the accretion 
rate of small dark matter halos of mass m into 
larger dark matter halos of mass M as a function 
redshift: dN / dt = A{z,M) F{m/M), where typi- 
cally A{z,M) oc (1 + z)2-2m0-15 and F{m/M) = 
{M/mY{l — m/MY . Fit parameters for merger 
rates in terms of dark halo mass, luminosity, or 
stellar mass are given in Table 1. 

2. The merger rate of galaxies of luminosities L > fL^, 
should evolve in a similar manner, with a redshift 
and luminosity dependence that follows /) oc 
(l-f/)(l + z)2-i. 

3. Only a small fraction (0.5% at z = 0, 10% at 
z = 4) of bright (> 0.4L*) galaxies should have 
experienced a major (> 0.3) merger in their very 
recent history (100 Myr, Figure 3 left panel). Even 
if mergers trigger the kind of short-lived, highly- 
efficient star formation bursts that are expected in 
some models, they cannot contribute significantly 
to the overall distribution of star formation rates 
at any given epoch. 

4. The predicted fraction of galaxies with a merger 
in the past 500 Myr, or alternatively within a past 
halo dynamical time, are in reasonable agreement 
with the fraction of galaxies that show observa- 
tional signs of morphological disturbance between 



Allgood, B., Flores, R. A., Primack, J. R., Kravtsov, A. V., 

Wechsler, R. H., Faltenbacher, A., & Bullock, J. S. 2006, 

MNRAS, 367, 1781 
Barnes, J. E. & Hernquist, L. 1996, ApJ, 471, 115 
Barton, E. J., Arnold, J. A., Zentner, A. R., Bullock, J. S., & 

Wechsler, R. H. 2007, ApJ, 671, 1538 
Bell, E. F., Phleps, S., SomerviUe, R. S., Wolf, C, Borch, A., & 

Meisenheimer, K. 2006, ApJ, 652, 270 
Bell et al. 2005, ApJ, 625, 23 

Berrier, J. C, Bullock, J. S., Barton, E. J., Gucnthor, H. D., 
Zentner, A. R., & Wechsler, R. H. 2006, ApJ, 652, 56 



redshifts z = — 2 (Figure 3, middle and left pan- 
els). We emphasize, however, that comparisons be- 
tween theory and observations suffer from signifi- 
cant uncertainties associated with mass-ratio de- 
pendencies and relaxation timescales. 

5. Galaxy merger rates should depend on at least 
three parameters: mass (or luminosity), merger 
mass ratio, and redshift (see Table 1). Therefore 
any attempt to compare two observational indica- 
tors of the merger rate or to relate specific obser- 
vations to theoretical predictions must take great 
care in the respective comparisons. 

6. Mergers must become increasingly important in 
shaping galaxy properties at z > 3. At z = 3, the 
fraction of galaxies with a merger in the past dy- 
namical time is ^ 4 times higher than at z = 0. We 
expect ~ 30% (60%) of > 0.4L* galaxies to have 
experienced a m/M > 0.3 major {m/M > 0.1 mi- 
nor) merger in the past 500 Myr at z = 3. Though 
it is unlikely that short-lived starbursts associated 
with these mergers drive the increase in the global 
star formation rate of galaxies with redshift, the 
broader implications of these mergers (fresh sup- 
ply of gas brought in to the central galaxy through 
accreted satellites, etc.) are undoubtedly linked to 
star formation and the general growth of galaxies 
on longer timescales. 
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